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Abstract. When a voltage is applied across a thin layer of cholesteric liquid crystal, fin- 
gers of cholesteric alignment can form and propagate in the layer. In computer simulation, 
based on experimental laboratory results, we demonstrate that these cholesteric fingers 
can solve selected problems of computational geometry, logic and arithmetics. We show 
that branching fingers approximate a planar Voronoi diagram, and non-branching fingers 
produce a convex subdivision of concave polygons. We also provide a detailed blue-print 
and simulation of a one-bit half-adder functioning on the principles of collision-based com- 
puting, where the implementation is via collision of liquid crystal fingers with obstacles 
and other fingers. 

1. Introduction 

Liquid crystals (LCs) are fluids composed of anisotropic (usually rod-shaped) molecules 
that exhibit long range order. In the simplest phase, known as the nematic, the molecules 
tend to align in a common direction, called the director n. In this paper we use a cholesteric 
phase, in which the LC contains at least one component which is chiral so that the director 
adopts a helical structure through the LC. In a display, a thin layer of LC is enclosed 
between transparent substrates which contain electrodes to enable the application of a 
voltage. 

LCs behave as elastic media and will adopt the director configuration that minimises the 
elastic energy of the system. In addition the molecules have anisotropic dielectric proper- 
ties so that applying an electric field will tend to cause them to rotate. The configuration 
adopted by the LC director will therefore be determined by the combination of the bound- 
ary conditions imposed by the inside surfaces of the substrates, as well as the interaction 
between the electric field and the LC molecules. In the devices considered in this paper 
the inside surfaces of the substrates are designed to induce homeotropic alignment, that is 
with the director orthogonal to the substrates. The thickness of the LC layer is chosen to 
be close to the natural pitch of the cholesteric LC, which is the distance over which the di- 
rector undergoes a full rotation. Under these conditions the natural tendency of the LC to 
twist is suppressed and the LC adopts a uniform untwisted configuration with the director 
orthogonal to the substrates. However, small perturbations to the system can upset this 
equilibrium and can cause the director to collapse into complex localised structures known 
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Figure 1. The surfaces are treated so that with no voltage apphed, the 
hquid crystal aligns with its director n aligned normal to the substrates, as 
shown in the left figure. When the voltage is applied, the electric field E 
rotates the LC, which has a negative dielectric anisotropy, so that it tends 
to align with the director parallel to the substrates. 



as cholesteric fingers. These fingers consist of extended domains of twisted LC and a rich 
variety of structures and phases has been reported [2T, ^26^, ^. 

In our system we perturb the equilibrium by applying a voltage. The LC that we use 
has a negative dielectric anisotropy, that is the dielectric constant is greatest for electric 
fields directed orthogonal to the long axes of the rod shaped molecules. The field, which 
is applied across the layer, tends to rotate the director away from the initial vertical state 
(Fig.[T]), and as the director becomes more planar the natural chirality of the LC dominates 
resulting in the formation of cholesteric fingers. The anisotropic optical properties of the 
LC means that these fingers are clearly visible under an optical microscope. In this paper 
we enhance the contrast by adding dichroic dyes [29j. The fingers nucleate at defects or 
particles [27] in the LC and grow as the applied voltage exceeds a threshold and retract as 
it reduces below the threshold. For larger voltages the fingers start to branch, and as they 
fill the space the fingers start to repel each other. 

Fingers have also been reported in other LC phases [19l |25] and there have also been 
studies of isolated fingers which can propagate as a result of electroconvection in nematic 
LCs [T5([23. These propagating localisations in liquid crystal first attracted the attention 
of the unconventional computation community in early 2000s in the frame of collision-based 
computing. It was demonstrated in cellular automata models that worms in liquid crystals 
can implement basic logical gates when they collide with each other and either annihilate 
or deflect as a result of the collision [2] however no extended results were obtained at that 
time. Recent work in HP Labs has focussed on the interaction between LC materials and 
engineered microstructures [301 EIJ some of these systems have shown the potential for 
engineering the nucleation and propagation of cholesteric fingers. These results inspired us 
to reconsider the concept of collision-based computing in liquid crystals and to expand a 
range of problems solved by the propagating and interacting fingers. 

The paper is structured as follows. In Sect. [2] we describe how we make our samples 
and measurements and in Sect. |3] we show how to imitate space-time dynamics of LC 
finger propagation using mobile automata on lattices. The classical problem of the planar 
Voronoi diagram is tackled by LC fingers in Sect. [Ij Subdivision of a concave shape in 
convex shapes using LC fingers is demonstrated in Sect. [Sj We describe a design of a one- 
bit half-adder implemented via collision between fingers and obstacles, and between pairs 
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of fingers, in Sect. [6| Final thoughts of future LC finger computing devices are outlined in 
Sect. [71 

2. Fabrications and measurement of devices 

Devices were constructed using glass substrates coated with a transparent conductor 
layer made from ITO. In order to control the nucleation of the cholesteric fingers we 
patterned polymer structures onto one surface. Fingers tend to nucleate at the sharpest 
points on any structures. Fig. [2] shows a typical set of features: 40/im diameter rings 
with 4 protrusions to seed the nucleation of fingers. The structures are S/xm high and 
cover the ITO on one substrate. They are made from SU8 (Microchem Corp.) using 
photolithography. The surfaces of the two substrates are then coated with a thin polymer 
layer to align the LC homeotropically. To assemble the cell, the substrates are gently placed 
in contact so that the SU8 structures set the cell spacing, and the substrates are sealed 
on two edges with UV curing glue (NOA73 from Norland Products Inc.). The cell is then 
capillary filled with the LC mixture, with the cell heated to above the isotropic/nematic 
phase transition temperature to avoid the flow artifacts that can occur when filling with the 
LC in the nematic phase. The LC used was MLC2037 (Merck KGaA) which was chosen as 
it has a negative dielectric anisotropy, and the filling was carried out at 95° C. It was doped 
with an additive (1.29% by weight of zliSll, Merck) to impart a chiral pitch to the mixture. 
The concentration of the chiral dopant was chosen so that the pitch was very close to the cell 
spacing, so that the homeotropic state was just stable. The mixture was then doped with a 
blend of dichroic dyes (1% by weight G232, G241 and G472 from Hayashibara Biochemical 
Laboratories, Inc.). These enhance the contrast and make the fingers visible, without the 
need to use polarisers [29j. Dichroic dyes are rod-shaped molecules that only absorb light 
polarised along their long axis. They align with the LC, so that in the quiescent state the 
LC and dyes are both aligned vertically. Thus the dyes absorb very little light. When the 
fingers form, the LC and dyes at least partially align at a more planar angle, so that the 
dyes absorb more light and the fingers appear black. Applying a voltage across the LC 
layer causes the fingers to nucleate from the sharp points of the structures. These then 
slowly extend while the voltage is maintained (Fig. [2]), and then retract when the voltage 
is removed. We use a IkHz sinusoidal voltage with an amplitude of around 1.5V. In the 
next section we describe a model that captures some of the behaviour of the cholesteric 
fingers. 

3. Mobile automata model of liquid crystal fingers 

A mobile automaton is a tuple A = {c, s, f,a,d,), where c G is an Euclidean co- 
ordinate of the automaton, c = {x,y); function / : R^ transforms coordinates as 
= -j-dsina and y^^^ = + dcosa; in models presented here d = 0.5. The automa- 
ton moves on a lattice L. All nodes of L are in state 'empty' (0) initially. At each step 
t of evolution time the nodes are updated as follows: for every u ^ L: = 1 if there is 
an automaton A with coordinates u. The 'occupied' state is absorbing. Obstacles are also 
defined as domains of 'occupied' states '1'. 



4 ADAMATZKY, KITSON, DE LACY COSTELLO, MATRANGA, YOUNGER 




(e) (f) 



Figure 2. Experimental laboratory snapshots (optical microscopy) of col- 
liding LC fingers, taken with an applied IkHz ac voltage with an amplitude 
of 1.5V. Each photo is taken at a different time: (a) taken before the voltage 
is applied, (b) 4s, (c) 25s, (d) 40s, (e) 55s and (f) 80s after the voltage is 
applied 
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Figure 3. Simulation of experimental results, shown in Fig. [2j using mobile automata. 
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Figure 4. Simulated interaction between fingers when distance be- 
tween initial seeds is large enough to allow for prolonged movement, 
(a) Solid/impassable boundaries, (b) Boundary-less space. 



Fig. [2] shows a sequence of microscope photos of fingers propagating away from a dense 
array of nucleation sites. As fingers approach each other they deviate from a straight path 
to avoid contact. We find that the fingers always tend to turn in the same direction, left 
in the case of figure 1. We believe this is due to the natural chirality of the LC mixture. 
The model is based on these observations. An automaton A chooses where to move as 
follows. If node {x^ + dsma,y^ ^dcosa) empty (0) then A moves into this node; otherwise 
it rotates left: a = a + 7r/360 and its state s is incremented. The exact angle of scattering 
is not known and may depend on many factors beyond our present knowledge. Thus we 
adopted increment 7r/360 which a is increased with until a 'free' site for the next position 
is found. The automaton stops when s exceeds a certain threshold of attempts. 

The mobile automata model is phenomenological. A mobile automaton moving in a 
two-dimensional discrete space very roughly imitates the tip of a LC finger propagating in 
a quasi-two-dimensional space being squashed between two electrodes. Our model in no 
way competes with existing numerical models of LC fingers, see e.g. [20j, but must rather 
be considered as a fast-prototyping tool in the design of LC finger computing algorithms. 

We qualitatively verified our model using experimental observations. For example if 
we consider Fig. [2j Discs with four singularities/protrusions are arranged in a regular 
grid so that the protusions are aligned between neighbouring discs (the protusions of four 
neighbouring discs describe a square array) (Fig. [2^). When a voltage of 1.5V is apphed, 
LC fingers nucleate near the protrusions (Fig.[2]3). The fingers propagate along the original 
axis determined by the protrusions. When two fingers come into a head-on collision, each of 
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them turns left (Fig.|2]3-f). To imitate this experimental finding we regularly and uniformly 
distributed clusters of mobile automata in a two-dimensional space (Fig. |3]). Initially the 
four automata of each cluster have their velocity vectors angles a oriented north, south, 
west and east (Fig. [3^). When finger f moving North collides with finger f" moving south 
the finger f turns south-east, and finger f" north-west (Fig. |3]3-d). Similarly, the finger 
initially travelling west deviates south-west and the finger travelling east turns north- 
east. Further changes of the fingers' trajectories are attributed to subsequent collisions 
(Fig. [3^f). If the distance between the initiation sites of the fingers is large enough to allow 
for prolonged movement of the fingers the four neighbouring fingers form spirals (Fig. ffl). 



4. Approximation of Voronoi diagram 

Let P be a non-empty finite set of planar points. A planar Voronoi diagram of the set P is 
a partition of the plane into such regions that, for any element of P, a region corresponding 
to a unique point p contains all those points of the plane which are closer to p than to any 
other node of P. A unique region vor{p) = {z G : d{p, z) < d{p, m) Mm G R^, m ^ z} 
assigned to the point p is called a Voronoi cell of the point p. The boundary of the Voronoi 
cell of the point p is built up of segments of bisectors separating pairs of geographically 
closest points of the given planar set P. A union of all boundaries of the Voronoi cells 
determines the planar Voronoi diagram: VD{P) — Up^pdvor{p) [23]. Voronoi diagrams are 
applied in many fields of science and engineering. A few books and conference proceedings 
are available on the theory and applications of the Voronoi diagram pTj [7], 

Construction of a Voronoi diagram is a classical problem of unconventional computing 
devices. This was the first ever problem solved in a reaction-diffusion chemical com- 
puter [21]. The basic concept of constructing Voronoi diagrams with reaction-diffusion 
systems is based on an intuitive technique for detecting the bisector points separating two 
given points of the set P. If we drop reagents at the two data points the diffusive waves, 
or phase waves if the computing substrate is active, travel outwards from the drops. The 
waves travel the same distance from the sites of origin before they meet one another. The 
points where the waves meet are the bisector points; see the extensive bibliography in [HH] 
and mechanisms of bisector formation in chemical media in [TOl [HI [121 [13 [HI • The Voronoi 
diagram is also approximated in crystallisation-based processors [5] and colonies of acellu- 
lar slime mould Physarum polycephalum [6J. Thus it is intuitive to ascertain how Voronoi 
diagrams could be approximated with LC fingers. 

The behaviour of the fingers is voltage-dependent. When the voltage exceeds a critical 
voltage (in this case 1.7 V)the fingers grow much more quickly and start to branch (Fig. [5]). 
Recursively branching fingers form fronts or a phase boundary between a free space and 
the domain filled by fingers. Fronts originating from different sources compete for space. 
When two or more fronts meet, they stop propagating. Thus to approximate a Voronoi 
diagram we place the seeding structures of fingers in positions of space corresponding 
to planar points from P and increase the voltage to cause finger branching. The loci of 
space not occupied by any fingers represents the edge of the Voronoi diagram VD(P). 
From experimental observation it can be concluded that the accuracy of Voronoi diagram 
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Figure 5. Experimental image of interaction between branching LC fingers 
originating from several seeds. The central point of the seeding structures 
are marked by small black discs. The edges of the Voronoi diagram are 
calculated using a classical sweepline algorithm |18j and are drawn as solid 
lines superimposed on the experimental image. 
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(d) (e) (f) 

Figure 6. Simulated approximation of Voronoi diagram, (a-e) Snapshots 
of the simulated dynamics of branching LC fingers. Edges of the Voronoi 
diagram calculated using a classical sweepline algorithm [ISJ are drawn by 
solid red/gray lines in (a-e) and shown explicitly in (f). 

construction is improved where fingers are highly branched (approximate a continuous 
expanding front emanating from a defined site). This can be achieved by controlling 
the voltage or maximising the density of fingers initiated by altering the geometry and 
positioning of the seeding structures (for example compare construction in Figs. [5^ and[5]3). 
The approximation of the Voronoi diagram of eleven planar sites in mobile automata model 
of LC fingers is shown in Fig. |6j 

5. Convex subdivision of concave polygons 

A Concave polygon is a shape comprised of straight lines with at least one indentation, 
or angle pointing inward. The problem is to subdivide the given concave shape into convex 
shapes. The problem is solved by LC fingers as follows. Fingers are generated at the 
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Figure 7. Convex subdivision of concave polygons in mobile automata 
model of LC fingers. Snapshots of the space after the propagation of all 
fingers has stopped, (a-d) Snapshots of LC finger model subdividing a 
concave polygon with three indentations, (e) Convex subdivision of star- 
shaped polygon. 

singular points of indentations and the fingers' propagation vectors are co-aligned with 
medians of the corresponding inward angles. Given a concave polygon every indentation 
initiates one propagating finger. A finger turns left, relative to its vector of propagation, 
when it collides with another finger or a segment of the polygon. Given a polygon with n 
indentations, n fingers will be generated. By following their "turn-left-if-there-is-no-place- 
to-go" routine and also competing for the available space with each other the fingers fill 
n — 1 convex domains. At least one convex domain will remain unfilled. 

In the example shown in Fig. [7^, there are three indentations. They initiate fingers 
propagating — approximately — south (a), west (6) and north-north-east (c). The finger 
travelling south (a) collides with a finger travelling north- north-east (c). In the result of 
the collision, finger a turn left and moves South-East. Finger c turns left and moves North- 
West- West (Fig.[7]b). At some stage finger b collides into finger c and finger a collides into 
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finger c. When fingers collide to the 'walls' of the polygon they turn left and more or less 
accurately follow the wall till next collision (Fig. [t]^). Each finger forms a spiral and gets 
stuck inside the spiral at some point of its growth. 

The computation of the convex subdivision of a concave polygon is completed when 
all fingers cease moving. The result of the computation is a set of convex polygons, each 
polygon is filled by a unique finger and at least one non-filled convex domain. In the 
example provided, a concave polygon is subdivided into four domains. The north-eastern 
domain is constructed by finger a, the south-eastern domain by finger c and the south- 
western domain by finger b. The north-western domain remains empty (Fig. [7]l). In the 
case of regular polygons, e.g. star-shaped in Fig. [7^, boundaries between finger-filled convex 
domains correspond to the skeleton of this planar polygonal shape. 



6. One-bit half-adder 

The half-adder presented is based on the paradigm of collision-based computing, which 
originates from conservative logic and billiard ball model by Fredkin and Toffoli [T7] and 
logical computation by colliding gliders in Conway's Game of Life [9j. A collision-based 
computer employs mobile compact finite patterns and travelling localisations to represent 
quanta of information in non-linear media. Information values, e.g. truth values of logical 
variables, are given by either absence or presence of the localizations or other parameters 
of the localizations. The localizations travel in space and when collisions occur the result 
can be interpreted as computation. There are no predetermined stationary signal channels 
(wires) — however stationary localisations, or reflectors, are allowed — a trajectory of the 
travelling pattern is considered a transient wire. Almost any part of the mediums' space 
can be used as a wire. Localizations can collide anywhere within the sample space, there 
are no fixed positions at which specific operations occur, nor location specified gates with 
fixed operations. The localizations undergo transformations, form bound states, annihilate 
or fuse when they interact with other mobile patterns. Information values of localizations 
are transformed as a result of colhsion [3]. 

In the design discussed we use reflectors. The reflectors are stationary structures or 
defects. They are 'artificial' in a sense that they do not belong to the liquid crystal medium 
but are externally and intentionally introduced to deflect propagating fingers. We assume 
the obstacles are rectangular. In real- world implementations corners of the rectangular 
reflectors would generate additional propagating fingers, thus obstacles must be smooth 
e.g. circular or oval. However, our model is coarse-grained and therefore any oval reflector 
would be composed of small rectangles and 'perceived' as a group of rectangles by the 
fingers. 

Basic gate xy is illustrated in Fig. |8| The gate consists of three obstacles oi, 02 and 03 
and two input fingers (Fig. [8^). The fingers represent values of Boolean variables x and 
y: presence of a finger x or y corresponds to x = 1 or y = 1 (True), absence of a finger 
corresponds to (False). Obstacles oi and 02 are used to deflect finger x while obstacle 
03 is for deflecting finger x after collision with finger y. 
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Figure 8. LC finger gate {x^y) {xy,xy). (a) Scheme: trajectories of 
fingers x and y entering gate are siiown by solid lines, trajectories of fingers 
after collision are shown by dotted lines, obstacles are shown by grey rect- 
angles and tagged oi, 02 and 03. (b-d) Configurations of two-dimensional 
mobile automata model. Fingers and obstacles are represented by states of 
cells, (b) Input x = 0, y = 1. (c) Input x = 1, ^ = 0. (d) Input x = 1, 
y = l. 



Let us consider input x = and y = 1 (Fig.jSp). Finger x is not present. Finger y travels 
south undisturbed. There are no obstacles in its way, thus it continues along its original 
trajectory which represents output xy. For input x — 1 and y — ^ only finger x enters 
the gate (Fig. [s]^). The finger x turns east after colhding with obstacle oi. It then collides 
with obstacles 03 and turns north. Both outputs are nil. When both inputs are True both 
fingers x and y enter the gate (Fig.jsjl). Finger x enters the gate early than finger y. Thus 
by the time finger y passes southward along obstacle 02 finger x moves northward. The 
fingers x and y collide with each other. In the result of this collision finger y turns east and 
finger x turns west. The deflected trajectory of finger y propagating east represents output 
xy. Deflected finger x colhdes with obstacle 02 deflected south, and eventually becomes 
trapped and propagation ceases (Fig. (sjl). 

To implement a binary half-adder one must realise calculation of a sum and carry on 
results: x ® y and xy. An architecture of LC finger half-adder can be implemented of 
three gates {x^y) {xy,xy) and five additional obstacles (Fig. [o]). There are seven input 
trajectories and two output trajectories. Inputs include three fingers representing constant 
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Figure 9. A scheme of one-bit half- adder implement able with LC fin- 
gers (a) and labels of key elements of the adder (b). 



Truth, two fingers representing values of x and two fingers representing values oiy. In the 
present model we do not tackle multiplication or splitting of signals therefore we assume 
copies of variables x and y are presented a priori. All input fingers travel south, output 
X ® y travels south and output xy travels east. Note that all input fingers enter the half- 
adder at different moments of time (Fig. [9^). In total we have seven fingers labelled a to 
g and fourteen obstacles labelled oi to 014 (Fig. (9)3). 

Simulation of one-bit half-adder in mobile automata model is shown is illustrated in 



Fig. 10 When both inputs x and y are False only three fingers corresponding to constant 



True propagate southward (Fig. MOa). Finger c collides is deflected by obstacle 04, turns 
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Figure 10. Configurations of two-dimensional mobile automata model. 
Fingers and obstacles are represented by states of cells, (a) x = 0, ^ = 0, 
(b) X = 0, y = 1, (c) X = 1, y = 0, (d) X = 1, ^ = 1. 
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east, is deflected by obstacle on and turns north. It collides with finger /. In the result of 
this collision finger c turns north-west and is "self-trapped" while finger / turns south-east. 
The body of finger / prevents further propagation of finger g. Finger / and g are deflected 
by obstacle 014 and propagate north. Thus for input x = and y = no output fingers 



For the input combination x = and y = 1 the situation develops as follows (Fig. [TOp). 

Fingers c, / and g^ representing constant True enter the the adder as usual. Fingers a and 
e, representing x, are absent. Fingers b and d, representing y = 1 enter the adder. Finger 
b collides with obstacle 03 and turns east, then collides with obstacle 07 and heads north. 
Thus fingers b and c come into a head-on collision. Finger b turns west and is "self-trapped" . 
Finger c turns east, collides with obstacle 07 and then oq and gets trapped. Finger d collides 
with obstacle os, heads east, collides with obstacle 012 and travels north. Fingers / and g 
continue their travel undisturbed. Thus operation (0, 1) (1,0) is imlpemented. 

When X — 1 and y — 1 (Fig. [To) :) finger a (representing first copy of x) is turned 
north by obstacles oi and 05. Finger c (representing second copy of x) is turned north by 
obstacles oio and 013, collides with finger /. Fingers c and / stop propagating in result of 
the collision. Finger c is turned east by obstacle 04 and then north by obstacle on. Finger 
g continues its propagation undisturbed and thus operation (1,0) (I5O). 

Most interactions between fingers take place in situation x — 1 and y — 1 (Fig. [To|d). 
Finger a collides with finger b. Finger b turns east as a result of the collision. Body of 
finger b prevent finger c propagating north. Finger d collides with finger e and gets "self- 
trapped". Finger e turns east in result of the collision and exits the adder. Final trajectory 
of finger e represents xy — 1 (Figs.|9^ and[To|d). Propagation finger g is blocked by body of 
fingers e and /, it does not exit the adder along its original trajectory, and thus x © y = 0. 



Cholesteric liquid crystals exhibit growth of localised phase defects — fingers — in a 
response to application of an ac electric field. The fingers show (almost) deterministic 
behaviour when they collide with other fingers and obstacles, thus they are suitable for 
implementation of collision-based computing schemes. For higher values of voltage applied 
fingers show branching. Wave-fronts of branching fingers stop their propagation when they 
collide with other fronts. Thus branching fingers can approximate Voronoi diagram, a 
plane subdivision based on proximity criteria. For low applied voltages the fingers remain 
solitary and thus each finger can represent a quantum of information and be an elementary 
unit of a collision-based computing device. To illustrate feasibility of the approach we 
provided design of a one-bit binary half-adder and proved correctness of its functioning 
using a mobile automata model. Collision-based computing prototypes presented in the 
paper are based on computer imitations of finger propagation and further work is required 
to implement the design in physical laboratory conditions. 





7. Discussion 



16 adamatzky, kitson, de lacy costello, matranga, younger 

References 

[1] Adamatzky A. Computing in Non-linear Media and Automata Collectives (loP Publishing, Bristol, 
2001). 

[2] Adamatzky A. New media for collision-based computing. In: Adamatzky A. (Ed.) Collision-Based 

Computing (Springer, 2003) 411-437. 
[3] Adamatzky A. (Ed.) Collision-Based Computing (Springer, 2003) 411-437. 

[4] Adamatzky A., De Lacy Costello B., and Asai T. Reaction-Diffusion Computers (Elsevier, Amsterdam, 
2005). 

[5] Adamatzky A. Hot ice computer. Phys. Lett. A 374 (2009) 264-271. 
[6] Adamatzky A. Physarum machines (World Scientific, 2010). 

[7] Anton F. In Proc. Sixth IEEE Int. Symp. Voronoi Diagrams, Copenhagen, Denmark, 2009. 

[8] Baudry J, Pirkl S, Oswald P. Looped finger transformation in frustrated cholesteric liquid crystals. 

Phys Rev E 59 (1999) 5562-5571. 
[9] Berlekamp E.R., Conway J.H., Guy R.L. Winning ways for your mathematical plays, vol. 2. Academic 

Press; 1982. 

[10] De Lacy Costello B. P. J. Constructive chemical processors — experimental evidence that shows that 

this class of programmable pattern forming reactions exist at the edge of a highly non-linear region. 

Int. J. Bifurcat. Chaos 13 (2003) 1561-1564. 
[11] De Lacy Costello B. and Adamatzky A. On multitasking in parallel chemical processors: experimental 

findings. Int. J. Bifurcat. Chaos 13 (2003) 521-533. 
[12] De Lacy Costello B. P. J., Hantz P., and Ratcliffe N. M. Voronoi diagrams generated by regressing 

edges of precipitation fronts. J. Chem. Phys. 120 (2004) 2413-2416. 
[13] De Lacy Costello B. P. J., Adamatzky A., Ratcliffe N. M., Zanin A., Purwins H. G., and Liehr A. The 

formation of Voronoi diagrams in chemical and physical systems: experimental findings and theoretical 

models. Int. J. Bifurcat. Chaos 14(7) (2004) 2187-2210. 
[14] De Lacy Costello B. P. J., Jahan I., Adamatzky A., and Ratcliffe N. M. Chemical tesselations. Int. J. 

Bifurcat. Chaos 19 (2009) 619-622. 
[15] Dennin M., Cannell D. S. and Ahlers G. Patterns of electroconvection in a nematic liquid crystal. 

Physical Review E 57 (1998) 638-649. 
[16] Dequidt A. and Oswald P. Does the electric Lehmann effect exist in cholesteric liquid crystals? Eur 

Phys J E 24 (2007) 157-166. 
[17] Fredkin E. and ToffoU T. Conservative logic. Int J Theor Physics 21 (1982) 219-253. 
[18] Fortune S. A sweepline algorithm for Voronoi diagrams. In Proc. 2nd Annu. Symp. Computational 

Geometry, Yorktown Heights, New York, 1986, pp. 313-322. 
[19] Li J.-F., Wang X.-Y., Kangas E., Taylor P. L, and Rosenblatt C, Suzuki Y., Cladis P. E. Reversible 

propagating fingers in an antiferroelectric liquid crystal. Phys. Rev. B 52, R13075R13078 (1995) 
[20] Nagaya T., Hikita Y., Orihara H. and Ishibashi Y. Numerical study of the growth of the cholesteric 

finger pattern J. Phys. Soc. Jpn. 65 (1996) 2713-2716. 
[21] Okabe A., Boots B., Sugihara K., and Chiu S. N. Spatial Tesselations (Wiley, Chichester, New York, 

2000). 

[22] Riecke H. and Granzow G. D. Localization of waves without bistability: Worms in nematic electrocon- 
vection. (1998) arXiv:patt-sol/9802003vl. 

[23] Preparata F. P. and Shamos M. I. Computational Geometry: An Introduction (Springer, New York, 
1985). 

[24] Tolmachiev D. and Adamatzky A. Chemical processor for computation of Voronoi diagram. Adv. 

Mater. Opt. Electron. 6 (1996) 191-196. 
[25] Wang X. Y., Li J.-F., Gurarie E., Fan S., Kyu T., Neubert M. E., Keast S. S., and Rosenblatt C. 

Kinetics of phase transition in an anticlinic liquid crystal induced by a iniform temperature field: 

Growth in one dimension Phys. Rev. Lett. 80 (1998) 4478. 



LIQUID CRYSTAL FINGERS COMPUTERS 



17 



[26] Pirkl S. and Oswald P. Travelling fingers in two-frequency cholesteric liquid crystals: fragments, loops 
and spirals Liquid Crystals 28 (2001) 299 

[27] Smalyukh LL, Senyuk, B.L, Palffy-Muhoray, P., Lavrentovich, O.D., Huang, H., Gartland, E.G., 
Bodnar, V.H., Kosa, T. and Taheri, B. Electric-field-induced nematic-cholesteric transition and three- 
dimensional director structures in homeotropic cells Phys. Rev. E 72 (2005) 061707 

[28] Ribiere, P., Oswald, P. and Pirkl, S. Crawling and spiraling of cholesteric fingers in electric field J. 
Phys. II France 4 (1994) 127 

[29] White, D.L. and Taylor, G.N. New absorptive mode reflective liquid-crystal display device J. Appl. 
Phys. 45 (1974) 4718 

[30] Kitson, S. and Geisow, A. Controllable alignment of nematic liquid crystals around microscopic posts: 

Stabilization of multiple states Appl. Phys. Lett. 80 (2002) 3635 
[31] Kitson, S.C., Edwards, E.G. and Geisow, A.D. Designing liquid crystal alignment surfaces Appl. Phys. 

Lett. 92 (2008) 073503 



